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A  model  for  the  oxidation  of  refractory  diborides 
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Wright-Patterson  AFB,  OH  45433-7817 
*  UES,  Inc.,  Dayton  OH  45432 
ANaval  Surface  Warfare  center,  Caderock,  MD 

Abstract 

Hypersonic  applications  require  materials  that  can  withstand  very  high 
temperatures.  Refractory  diborides  appear  to  be  the  most  promising  candidates  for 
these  applications,  from  the  standpoint  of  resistance  to  oxidation  and  evaporative  loss. 
In  this  work,  a  mechanistic  model  is  presented  that  predicts  the  oxidation  behavior  of 
diborides  of  Zr,  Hf  and  Ti.  Using  available  thermodynamic  data  and  literature  data  on 
vapor  pressures,  good  correspondence  is  obtained  between  theory  and  experiments  for 
weight  gain,  recession  and  scale  thickness  as  a  function  of  temperature  and  oxygen 
partial  pressure. 

1.0  Introduction 

The  leading  edge  of  a  hypersonic  vehicle  is  subjected  to  very  high  temperatures 
and  has  to  survive  the  flight  without  losing  the  integrity  of  the  leading  edge  geometry 
along  the  span.1,2  The  high  temperatures  and  the  low  pressure  environment  make 
oxidation  and  evaporative  loss  critical  factors  in  material  choice.3  The  diborides  of 
refractory  metals  (esp.,  Zr  and  Hf),  which  are  termed  ultra-high  temperature  ceramics 
(UHTC)  have  high  melting  points  and  are  the  most  promising  among  known  materials 
that  are  likely  to  survive  for  long-term  use  under  hypersonic  conditions.3,4  The  diborides 
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are  also  attractive  due  to  their  higher  thermal  conductivity  which  will  be  required  to  offset 
the  stringent  thermal  gradients  imposed  by  the  aerothermal  effects  of  hypersonic 
conditions. 

The  history  of  work  on  oxidation  of  diborides  is  given  by  Opeka  et  al.3  After  the 
promise  shown  from  early  works,  a  good  screening  study  of  refractory  element  diborides 
was  conducted  and  reported  by  Fenter.5  He  reported  the  diborides  of  Zr  and  Hf  to  be 
the  most  oxidation  resistant  of  the  diborides.  He  further  reported  that  additions  of  SiC  up 
to  20%  by  volume  further  improved  the  oxidation  resistance  of  these  materials.  It  has 
been  pointed  out  by  Opeka  et  al. 3,4,6  that  the  high  melting  point  and  low  vapor  pressures 
of  the  oxides  and  sub-oxides  of  Zr  and  Hf  are  responsible  for  the  superior  high 
temperature  resistance  of  these  materials.  Using  volatility  diagrams,  they  further  point 
out  that  of  the  various  possible  high  melting  compounds  of  Zr,  and  Hf,  the  borides 
appear  to  have  the  least  issues  with  respect  to  disruption  of  the  scale  formation.  In 
carbides,  CO  vapor  pressures  exceed  1  atm  above  1730C,  and  in  silicides,  the  active 
oxidation  to  gaseous  SiO  disrupts  the  scale  above  1800C.  In  contrast,  the  boric  oxide  is 
glassy,  flows  to  fill  porous  oxide  scale,  allowing  evaporation  of  the  oxide  to  take  place 
without  breaking  up  the  scale.  When  SiC  is  added  to  the  system,  further  improvements 
are  possible  since  the  silica  based  glass  offers  greater  resistance  to  evaporation  than 
B203  alone.7 

Thus  far  understanding  of  the  superior  oxidation  behavior  of  the  diborides  has 
been  qualitative.  For  engineering  design  and  application,  a  quantitative  model  that 
predicts  all  aspects  of  oxidation  is  desired.  It  is  desirable  to  have  predictions  under  a  set 
of  complex  conditions  of  temperature  and  environment  of  factors  such  as  scale 
thickness,  metal  recession  and  weight  change.  Many  experimental  work  have  been 
under  isothermal  conditions  and  in  laboratory  environment,  while  real  engineering 
conditions  will  involve  thermal  gradients  and  environmental  gradients. 
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A  model  that  can  be  applied  to  engineering  conditions  is  the  long-term  objective 
of  this  work.  However  the  complex  nature  of  the  materials  make  it  necessary  to  start 
with  a  model  for  the  simpler  systems  and  conditions,  where  reliable  experimental  data 
under  well  known  conditions  are  available.  In  this  paper,  we  present  a  model  that  is 
found  to  predict  the  scale  thickness,  metal  recession  and  weight  changes  for  diborides  of 
Zr,  Hf,  and  Ti  under  isothermal  conditions.  Various  factors  such  as  residual  boria 
content  and  dependences  on  temperature  and  oxygen  partial  pressures  are  also 
predicted  and  found  to  be  in  reasonable  agreement  with  experiments.  It  is  suggested 
that  this  model  provides  a  reasonable  basis  for  future  work. 

2.0  Conceptual  framework 

The  model  reported  in  this  work  is  based  on  concepts  of  the  oxidation  process 
derived  from  microstructural  details  available  in  the  literature.  The  morphology  of  the 
oxidation  product  as  well  as  the  phase  present  and  their  distribution  have  been  reported 
on  by  several  works.  For  the  case  of  ZrB2  (and  HfB2),  oxidation  in  air  at  temperatures  of 
up  to  1300°C  results  in  a  dense  adherent  oxide  scale.  The  scale  consists  of  a  Zr02  and 
B203,  and  no  intermediate  phases  have  been  observed  at  the  interfaces.  The  top 
surface  of  the  scale  is  nodular  Zr02,  and  the  cross-section  reveals  more  columnar 
structure  of  Zr02  with  glassy  B203  in-between  the  columnar  grains.  At  lower 
temperatures  (<1 100C)  a  glassy  B203  is  observed  on  top  of  the  Zr02  scale,  but  is  absent 
above  1300C.  The  microstructure  of  the  Zr02  appears  to  change  from  equiaxed  grains 
at  temperatures  below  1 10C,  to  columnar  at  higher  temperatures.  This  transition 
temperature  is  ~1650C  for  HfB2.  Similar  transition  has  been  noted  with  respect  to 
oxygen  partial  pressure.  The  oxygen  partial  pressure  dependence  changes  from  linear 
to  no  dependence  as  temperature  is  increased  above  1 150C. 
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Based  on  the  above  observations,  the  following  conceptual  framework  is 
suggested.  A  schematic  sketch  is  shown  in  Fig.  1 .  Oxidation  of  the  refractory  element 
results  in  an  evaporation  resistant  refractory  oxide  which  forms  a  porous  skeleton.  The 
oxidation  of  boron  results  in  a  glassy  boric  oxide  that  flows  to  fill  the  porous  skeleton.  At 
the  surface  the  boria  evaporates.  At  steady  state  the  oxidation  takes  place  through 
gaseous  diffusion  of  oxygen  through  the  pores  in  the  oxide  skeleton,  continuing  to 
diffuse  through  the  liquid  boria  (that  fills  the  porous  skeleton  beneath  the  surface)  to 
reach  the  ZrB2  oxide  scale  interface.  The  two  rate  limiting  reactions  are 

ZrB2  +  “O2 —  >  Zr02  +  B202(l )  ...  (1 ) 

B203  (I)  ->  B203  (g)  ...(2) 

The  rate  of  dissolution  of  oxygen  in  B203  is  taken  not  to  be  rate-limiting  in  the 
model  developed  below.  The  permeation  of  oxygen  through  the  refractory  oxide  is 
assumed  to  be  negligible  compared  to  that  through  the  pores  filled  with  glassy  boria. 

3.0  The  Model 

The  objective  of  the  model  is  to  predict  the  kinetics  of  scale  thickness,  substrate 
recession  and  weight  change  as  a  function  of  time  at  temperature  and  oxygen  partial 
pressure.  A  schematic  of  the  model  is  shown  in  Fig. 2. 

The  ambient  air  is  taken  to  consist  of  N2  and  02  only.  A  small  fraction,  fc,  of  the 
scale  is  taken  to  be  porous  and  provide  a  continuous  pathway  for  liquid  B203  to  pass 
through  and  evaporate  at  the  surface,  which  is  taken  to  be  a  perfect  sink  for  B203.  The 
relative  kinetics  of  filling  of  the  pore  by  oxidation  at  the  scale/substrate  interface  and  the 
evaporation  at  the  surface  of  B203  determines  the  extent  to  which  the  pores  are  filled 
with  B203.  At  a  given  time,  t,  the  scale  has  a  thickness  of  L  and  the  inner  portion  of  the 
scale  is  filled  with  liquid  B203  to  a  thickness,  h.  The  distance  (L-h)  is  the  diffusion 
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distance  over  which  gaseous  oxygen  must  diffuse  to  get  to  the  B203.  The  oxygen  is 
taken  to  dissolve  in  B203  and  continue  to  diffuse  to  the  substrate/scale  interface.  The 
oxide  scale  itself  is  taken  to  be  impervious  to  oxygen  (due  to  its  low  electrical 
conductivity). 

At  the  scale/substrate  interface,  reaction  1  takes  place.  Thus,  at  steady  state, 
the  flux  of  oxygen  must  be  balanced  by  the  formation  rates  of  Zr02  and  B203. 


5  5 

7  =  —  7  =  —  7 

J  02  2  S2°3  2  Zr°2 


...  (3) 


The  flux  of  oxygen  down  the  porous  channel  that  is  (L-t)  in  length  is  given  by  : 


Jn  =  Dn 


Ca  -C 
*-"0 

L-h 


■fc 


...  (4) 


In  the  above  equation,  D  refers  to  diffusion  coefficient,  C  refers  to  concentration, 
subscript  02  refers  to  oxygen,  and  superscript  "i"  refers  to  the  interface  between  B2O3  (1) 
and  B203(v),  which  is  at  a  distance  t  from  the  substrate/scale  interface.  Similarly  the  flux 
ofB203(v)  from  the  surface  of  B203(1)  to  the  ambient,  is  given  by  : 


•fso  —  D 


C 


B?0-, 


C 


BjO-i 


B7Oi 


L-h 


fc 


(5) 


Combining  equations  (3),  (4)  and  (5)  the  oxygen  concentration  at  the  B203  liquid  vapor 
interface  is  obtained  as  : 


>0/  _  a 

U o7  ~  ~ 


5  A 


2  Dn 


(Qo,  -CU) 


(6) 


The  concentration,  C^0} ,  at  the  ambient  surface  is  taken  to  be  zero.  The  concentration, 
,  of  the  B203(v)  at  this  interface  is  obtained  from  vapor  pressure  data.  8 


C' 

^ B20 3 


lB2o} 

RT 


PBO  =  3x1 010  exp(-^’^^)  atm  ...  (7) 
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The  permeability  of  oxygen  through  the  liquid  B2O3  depends  on  the  oxygen 
partial  pressures  across  the  liquid  layer  of  thickness  t.  The  oxygen  partial  pressure  at  the 
B2O3  liquid-vapor  interface  is  given  by, 

Po,=RTC'0i  ...(8) 


The  partial  pressure  of  oxygen  at  the  scale  B203  interface  is  given  by, 

a  P5/2 

ts  _  ]_ZrB2r_  02  /Q\ 

^  reaction 1  *  *  *  \^ ) 

aB203aZr02 

O 

Using  the  thermodynamic  data  of  Barin  for  reaction  1,  equation  (9)  gives  the  following 
for  the  oxygen  partial  pressure,  P£ ,  at  the  ZrB2-B203  interface  (i2). 


Pq2  =5x1  O10 exp( - — — )  atm  ...  (10) 

The  values  given  by  equation  (10)  are  consistent  with  the  maps  presented  by 
Fahrenholtz.9  Diffusivity  of  oxygen  through  boria  has  been  measured  by  Tokuda  et  al., 
who  found  a  liner  dependence  on  oxygen  partial  pressure  and  deduced  that  it  diffuse  as 
molecular  oxygen.  The  oxygen  flux,  J0i(B20 3)  >  across  the  liquid  boria  is  related  to  the 

partial  pressure  gradient  through  oxygen  permeability  coefficient,10’11  as: 

j  ^02ZB203  rr pi  pi 2i  Q  If 

t/02(B203)  ^  J  cl1  02  1  02  J  •"  V11/ 

The  oxygen  permeability  coefficient,  P01_Bim ,  is  obtained  from  the  literature11’12,  as: 


02-B203 


n  ,n-3  14000  ,  . 

=  7x10  exp( - )mole/ m-  s 


■  atm 


(12) 


Combining  equations  1 1  and  12  and  equating  to  the  flux  of  oxygen  through  the  gaseous 
layer,  given  by  equation  4,  an  equation  relating  the  boria  layer  thickness,  t,  to  the  scale 
thickness,  L,  is  obtained. 
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h  =  qL 


_ Pq2-B2Q3  C^Q2  P()2  ) _ 

*>0,  (Co,  ~  Cq,  )  +  Poi-BlOiC (Pq2  ) 


...(13) 


The  equations  for  the  rate  of  change  of  scale  thickness  and  mass  of  the  scale  are 
given  by  accounting  for  mass  balance  and  equation  3. 


dL  1  dW7r(h  1 


12  ,  M7 


dt  1  -fc  dt  pZrQ  l-/c  5  2  pZr„ 


Combining  with  equation  4,  one  obtains: 


dL  _  2  M  zro,  d  Co2  Cq2  /' 
dt  5  pZr0  °2  L(l-q )  1  -  fc 


...  (15) 


Integration  of  equation  15  gives  the  parabolic  growth  equation  for  scale  thickness. 


L  =  2[— 


2  MZro2  Cq2  Cq2  fc 


(i -q)  i -fe 


...  (16) 


The  recession  of  the  substrate  is  given  by 


R=L(l~fc) 


M ZrB2  !  PlrBn 
M  ZrO,  I  PzrO, 


...  (17) 


The  total  weight  change  per  unit  area  is  given  by: 


—  =  LpZrO,  (1  -fc)  +  ¥cPb203  -  RpZrB2 


...  (18) 


It  remains  to  obtain  expressions  for  the  diffusion  coefficient  of  oxygen  through  the 
porous  region  filled  with  oxygen,  ECCflv)  and  nitrogen.  The  diffusion  coefficient  in  a 
multicomponent  gaseous  system  can  be  approximated  by:13 


'1,(2, ..i)  M 


Yj(Xi'D 1-1 ) 


X<  =  77T 


...  (19) 
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Where  D  refers  to  diffusivity,  n  to  mole  fraction,  subscript  "1,  (2„i)"  refers  to  the 
diffusivity  of  specie  1  in  a  medium  of  i  species,  subscript  "1-i"  refers  to  diffusivity  of 
specie  1  in  a  binary  mixture  of  species  1  and  i.  The  diffusivity,  Di_2,  of  specie  1  in  a 
binary  gas  mixture  with  specie  2,  is  given  by14 


0.0018583T3/2  Vd/Ml)  +  (1/M  2) 
Pr2  nD 


...  (20) 


where  D]  _2  =  gas  diffusivity  in  cm2  /  sec 

_  1.06036  (  0.193000  |  1.03587  |  1.76474 

D  ~  T*015610  +  exp(0.47635T*)  +  exp(1.52996T*)  +  exp(3.8941  IT*) 
kT  I - 

T*  =  — ,  sn  =  ,  ru  =  0.5(rj  +  r2)  in  Angstroms,  P  =  pressure  in  atm, 

£n 

Mt  =  molecular  weight  (gms  /  mole) 


The  parameters  needed  for  the  above  expressions  can  be  obtained  from  the  work  of 
Svehla.15  When  the  size  of  the  pore  is  small,  the  diffusivity  of  the  gases  is  determined  by 
Knudsen  diffusivity.  This  is  given  by16: 


_4(8RT 

4  V 


7lM 


r 

2 


..(21) 


where  Dk  is  the  Knudsen  diffusivity,  M  the  molecular  mass  of  the  diffusing  species,  and  r 
the  radius  of  the  porous  pathway.  The  effective  diffusivity  is  given  by16: 

Deff  =(AT‘  +A-2"1)"1  -  (22) 


Using  the  expressions,  19  through  22,  Equations  16  through  18  give  the  scale 
thickness,  recession  and  weight  gain  as  a  function  of  temperature,  time  and  partial 
pressure  of  oxygen. 
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4.0  Model  Predictions  compared  with  literature  data 

4.1  ZrB  7 

The  model  was  verified  by  comparing  with  experimental  data  reported  in  the 
literature  for  ZrB2.  Extensive  TGA  data  on  weight  change  versus  time  in  250Torr  of  pure 
oxygen  has  been  reported  by  Tripp  and  Graham.17  They  also  compare  their  results  with 
the  parabolic  rate  constant  data  from  prior  work  by  Berkowitz  and  Kuraikose  and 
Margrave.19  Figure  3  shows  their  results  compared  with  the  model  predictions,  for  an 
assumed  pore  fraction  of  0.05  and  pore  radius  of  0.5  micron.  The  model  predictions  for 
the  weight  change  with  time  are  shown  for  three  different  temperatures  in  Figure  3(a), 
while  the  parabolic  rate  constants  are  shown  plotted  in  Fig3(b).  The  data  of  Tripp  and 
Graham  deviates  from  the  model  at  temperatures  above  1200C.  However  when  the  model 
is  compared  with  all  of  the  reported  data,  the  correspondence  is  good. 

A  study  by  Fenter  reported  on  the  metal  recession  of  several  diborides  in  air,  as  a 
function  of  temperature.5  The  results  on  ZrB2  is  shown  in  Figure  4(a)  along  with  the 
model  predictions.  The  model  compares  very  well  with  data  upto  -1850C.  Above  this 
temperature,  the  model  predicts  that  all  of  the  boria  will  evaporate  as  soon  as  it  forms. 

The  experimental  data  shows  a  significant  enhancement  in  recession  at  this  temperature. 
The  model  predicts  the  transition  temperature  correctly,  but  not  the  upward  trend  in 
recession.  It  is  not  clear  if  there  was  significant  spallation  in  the  experiments  above  this 
temperature. 

Finally,  recent  works  by  Opeka  et  al. 3,6  and  Talmy  et  al.20  have  reported  on  mass 
change  and  oxide  layer  thicknesses  of  ZrB2  in  air  or  Ar/C>2  mixtures.  Their  data  are 
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shown  compared  with  the  model  predictions  in  Figure  4  (b)  and  (c).  The  correspondence 
is  seen  to  be  good.  The  choice  of  pore  fraction  and  radius  were  kept  the  same  for  all  of 
the  model  predictions. 

4,2  HfB? 

The  model  is  easily  extended  to  other  diborides  by  simply  using  the  appropriate 
thermodynamic  data  and  the  appropriate  physical  properties  such  as  molecular  weight, 
density,  etc.  For  HfB2,  one  obtains  the  oxygen  partial  pressure,  P^ ,  at  the  HIB2-B2O3 
interface  as 

P'0](m  =  6x1 09  exp(-  1QQ^>7Q)  atm  ...  (23) 


There  have  been  two  reports  in  the  literature  that  have  characterized  the  oxidation 
behavior  of  FHB2.  Fenter5  has  measured  the  recession  rates  in  air,  while  Berkowitz18  has 
reported  on  the  parabolic  rate  constants  in  250Torr  of  pure  oxygen.  These  data  are 
compared  with  the  model  predictions  using  a  pore  fraction  of  0.05  and  pore  radius  of  0.5 
micron  in  Figure  5.  The  predicted  recession  is  higher  than  the  measured  values,  while  the 
predicted  parabolic  rate  constants  are  lower. 

4,3  TiB? 

For  TiB2,  one  obtains  the  oxygen  partial  pressure,  Pq  ,  at  the  TiB2-B203  interface 
as 

P‘l{Hf  )  =  5xl010  exp(-^^)  atm  ...  (24) 
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Using  the  appropriate  physical  properties  for  Ti02  and  TiB2  the  oxidation  behavior  was 
predicted,  once  again  assuming  a  pore  fraction  of  0.05  and  a  pore  radius  of  0.5  micron. 
The  predictions  were  evaluated  by  comparing  with  experimental  data  available  on  weight 
change  and  recession  as  reported  by  Koh  et  al.21  This  comparison  is  shown  in  Figure  6. 

5.0  Parametric/sensitivity  studies 

A  parametric  study  was  conducted  to  assess  model  sensitivity  and  to  determine 
the  effects  of  experimental  variables.  Figure  7  shows  the  temperature  dependence  of 
the  parabolic  rate  constant  in  ZrB2.  A  clear  change  in  apparent  activation  energy  is  seen 
at  around  1400C.  The  figure  includes  the  prediction  fraction  of  boria  in  the  scale.  It  is 
clear  that  the  change  in  activation  energy  arises  from  the  loss  of  boria  from  the  scale. 

Next,  the  model  was  used  to  study  the  effect  of  varying  the  pore  fraction  and 
pore  radius,  to  assess  its  sensitivity.  These  parameters  are  difficult  to  measure  and 
have  not  been  reported  experimentally.  Further,  since  the  model  assumes  neglects  the 
tortuosity  of  the  porous  scale,  the  pore  radius  and  and  pore  fraction  used  in  the  model 
must  be  considered  to  be  "effective"  parameters  rather  than  the  actual  physical  values. 

It  is  quite  possible  that  these  parameters  themselves  vary  with  temperature  or  other 
experimental  conditions  such  as  environment.  Thus  it  is  useful  to  conduct  a  sensitivity 
study. 

Figure  8  shows  the  effect  of  pore  radius  on  the  recession,  scale  thickness, 
parabolic  rate  constant  and  the  fraction  of  boria  in  the  scale.  It  is  found  that  the  effect  is 
significant  when  the  pore  radius  is  smaller  than  about  1  micron.  Further,  the  effect  is 
insignificant  at  temperatures  below  about  1400C,  which  corresponds  to  the  temperature 
where  boria  evaporation  becomes  significant. 

Figure  9  shows  the  effect  of  pore  fraction  on  the  oxidation  behavior  of  ZrB2  in  air. 
The  rate  constant,  scale  thickness  and  recession  all  increase  by  about  an  order  of 
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magnitude  when  the  pore  fraction  is  increased  from  0.025  to  0.2.  Thus  the  effect  is 
nearly  linear. 

Figure  10a  summarizes  the  effect  of  oxygen  partial  pressure  on  the  parabolic 
rate  constant  as  a  function  of  temperature.  The  dependence  is  linear  at  1000C  and 
becomes  independent  at  1600C;  this  is  consistent  with  the  experimental  observation  of 
Berkowiz.18  Figure  10b  shows  that  pore  radius  is  predicted  to  have  a  significant  effect  at 
low  oxygen  partial  pressures  and  the  pore  fraction  has  a  near-linear  effect  as  might  be 
expected. 

6.0  Summary 

A  physical  model  to  predict  the  oxidation  behavior  of  refractory  diborides  is 
presented.  The  model  assumes  that  the  refractory  oxide  is  inert  to  oxygen  diffusion  due 
to  lack  of  sufficient  electrical  conductivity.  The  porosity  in  the  oxide  is  taken  to  be  filled 
with  liquid  boria  which  evaporates  from  the  surface.  A  comparison  with  experimental 
data  in  the  literature  shows  that  the  steady  state  model  is  found  to  predict  the  weight 
change,  metal  recession,  oxide  scale  thickness,  the  temperature  dependence  of  the 
parabolic  rate  constant  reasonably  well  for  ZrB2.  Comparison  with  limited  literature  data 
on  HfB2  and  TiB2  show  that  the  model  may  have  a  general  applicability  to  refractory 
diborides,  although  data  on  electrical  conductivity  of  the  oxides  especially  at  high 
temperatures  will  be  required  to  refine  the  model. 
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Figure  1  Schematic  sketch  of  mechanisms  involved  in  the  oxidation  of  ZrB2  in  air,  that 
were  assumed  in  the  model. 
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B203(l)  formation  rate  =  Flux  of  B203(g)  OUT  =  (  5/2)  x  Flux  of  02  (g)  IN 


02  permeation  through  B203  (I)  =  Flux  of  02  (g)  IN 


air 

(0.2  02,  0.8  N2) 

perfect  sink 
for  B203(g) 


PB2o3  (at  l-g  interface) 

=  Vapor  Pressure  in  eqbm  with  B203(l) 


Figure  2  A  schematic  sketch  of  the  steps  involved  in  the  oxidation  of  ZrB2,  that  were 
considered  in  the  model.  The  Zr02  is  taken  to  be  inert  to  oxygen  permeation. 
The  porous  region  within  the  Zr02  is  assumed  to  be  continuous  and  tortuosity  is 
neglected.  The  pores  were  taken  to  be  filled  with  the  oxidation  product,  liquid 
boria,  which  is  lost  partially  by  evaporation  through  the  porous  region.  At  steady 
state  the  porous  region  is  partially  filled  with  boria. 
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1  /  Temperature,  K 

1  7 

Figure  3  (a)  The  experimental  data  from  Tripp  and  Graham  ,  on  the  measured  weight 
gain  of  samples  in  250  Torr  of  pure  oxygen  at  different  temperatures  compared 
with  the  model  predictions,  (b)  The  predicted  parabolic  constant  as  a  function  of 

1  7 

temperature  is  compared  with  the  data  of  three  different  works  in  the  literature. 

19  The  porosity  in  the  scale,  fc,  was  taken  to  be  0.05  and  the  pore  radius  to  be  0.5 
micron. 
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Figure  4  (a)  Model  predictions  of  the  recession  in  air  as  a  function  of  temperature  is 

shown  compared  with  data  of  Fenter,  1 97 1 ,5  (b)  The  model  predictions  for  mass 
change  is  shown  compared  with  experimental  data  obtained  in  Ar/02  gas  in  2 
hrs.6,20  (c)  The  scale  thickness  obtained  after  5hrs  in  air  by  Opeka  et  al.6  is  shown 
compared  to  the  predictions.  The  porosity  in  the  scale,  fc,  was  taken  to  be  0.05 
and  the  pore  radius  to  be  0.5  micron. 
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Figure  5  Comparison  of  model  predictions  with  experimental  data  on  HfB2.18  (a) 

recession  in  flowing  air  from  Fenter,  1971  and  (b)  Parabolic  constant  for  weight 
gain  as  a  function  of  temperature  from  Berkowitz,  1966,  in  25  OT  of  oxygen.  The 
porosity  in  the  scale,  fc,  was  taken  to  be  0.05,  and  pore  radius  to  be  0.5  micron. 
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Figure  6  Comparison  of  model  predictions  with  experimental  data  on  TiB2.  (a)  weight 

gain  with  time  in  air  and  (b)  scale  thickness  after  oxidation  in  air.  The  porosity  in 
the  scale,  fc,  was  taken  to  be  0.05,  with  a  pore  radius  of  0.5  micron.  The  data 
were  obtained  from  Koh  et  al,  2004.21 
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Figure  7  The  model  predicts  a  change  in  the  apparent  activation  energy  of  the  parabolic 
rate  constant  that  correlates  with  loss  of  boria  from  the  scale. 
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Figure  8  Sensitivity  studies  show  that  the  pore  radius  (Knudsen  effect)  has  a  significant 
effect  when  it  is  below  about  1  micron.  The  effect  on  (a)  recession,  (b)  boria 
fraction  in  scale,  (c)  oxide  thickness  and  (d)  parabolic  rate  constant  are  shown. 
The  effects  are  insignificant  below  1400C,  where  oxygen  diffusion  through  boria 
is  rate  limiting. 
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Figure  9  Sensitivity  studies  show  that  the  effect  of  pore  fraction  is  nearly  linear.  It  is 
seen  that  there  is  an  order  of  magnitude  change  in  rate  constant,  oxide  thickness 
and  recession,  for  a  change  in  pore  fraction  from  0.025  to  0.2. 
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Figure  10  (a)  Parametric  studies  show  that  the  dependence  of  the  parabolic  rate  constant 
on  the  oxygen  partial  pressure,  changes  from  unity  at  1000C  to  zero  above  1600C. 
The  P02  dependence  itself  is  dependent  on  (a)  the  pore  size  and  (b)  pore  fraction. 
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